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摘 要 为 本 文 基于 热 质 耦 合 的 Lattice Bhatnagar-Gross-Krook(CLBGK) 模型 ， 通 过 引入 浓度 分 布 函 数 ， 利 用 格子 
Boltzmann 方法 对 顶 盖 驱 动 的 复合 方 腔 内 的 双 扩散 混合 对 流 现象 进行 了 研究 ,复合 方 腔 由 多 和 孔 介 质 区 域 和 自由 空间 组 成 
分 析 了 路 易 斯 数 Le=2.0, 浮 升 力 比 N=1.0, 格拉 晓 夫 数 Gr=104 和 普 朗 特 数 Pr=0.7 时 , FLBRA (e = 0.6/0.7/0.8)、 方 
腔 中 多 孔 介质 层 位 置 及 理 查 德 森 数 Ri(107? <Ri<103) 对 内 部 混合 对 流 及 热 质 扩散 的 影响 . 给 出 了 方 腔 内 温度 、 浓 度 和 
流 线 分 布 ， 以 及 高 温 高 浓度 壁面 的 平均 努 塞 尔 数 Nuav 和 平均 舍 伍 德 数 Shav。 研究 结果 表明 : 多 和 孔 介 质 层 对 顶 盖 驱 动 方 
腔 内 热 质 双 扩散 影响 显著 ， 且 方 腔 左 壁 壁面 平均 努 塞 尔 数 Nua SPIRER Shav， 在 位 置 Di~D3 之 间 随 多 孔 介质 
层 的 右 移 而 增 大 , 在 位 置 D3 上 随 Ri(10 <Ri<103) 的 增 大 而 减 小 . 


对 流 格子 
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Lattice Boltzmann Simulation of Double Diffusive Mixed Convection in a 
Lid-Driven Composite Enclosure 
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(School of Energy and Power Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China) 


Abstract Based on the coupled lattice Bhatnagar-Gross-Krook (CLBGK) model, this article adopts 
lattice Boltzmann method and introduces concentration distribution function to investigate the double 
diffusive mixed convection in a lid-driven composite enclosure, which is composed of the space filled 
with porous medium and free space of pure fluid. At Le=2.0, N=1.0, Gr=10*, Pr=0.7, the influence 
of porosity (e = 0.6/0.7/0.8), porous medium layer locations D, and Richardson number Ri(10~° < 
Ri<103) on the double diffusive heat and mass transfer is performed. Results are presented in terms 
of isotherms, isoconcentrations, streamlines and the average Nusselt and Sherwood numbers along 
the left high temperature and concentration wall. Results show that the porosity and locations of 
porous medium have a significant influence on the doublediffusiveheat and mass in the lid-driven 
enclosure. The average Nusselt and Sherwood numbers along the left wall increase with the mov- 


ing right of porous medium layer from D1 to D3 and decrease with the increasing of Ri from 107 to 10°. 


Key words composite enclosure; lid-driven; double diffusion; lattice boltzmann simulation 
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多 孔 介 质 与 纯 流 体 空间 组 合成 的 复合 区 域 中 的 。 进行 了 广泛 研究 . 如 流体 为 达 西 流 或 非 达 西 流 时 , 多 


质 双 扩 散 现 象 是 自然 界 和 工农 业 生 产 中 普遍 存在 
的 一 种 现象 ， 如 室内 有 机 挥发 性 化 合 物 (VOCs) 的 
扩散 问题 、 土 壤 中 的 污染 传播 问题 、 农 业 生 产 中 的 
干燥 问题 ， 化 工 食品 加 工 中 的 脱水 及 晶体 生长 问题 
等 。 因此， 对 复合 区 域内 的 热 质 双 扩 散 现象 的 研究 
具有 十 分 重要 的 意义 。 

学 者 们 对 复合 方 腔 内 的 传 热 与 热 质 双 扩散 现象 
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孔 介 质 层 的 物性 参数 (如 : LRE, HER, LIT 
质 层 的 厚度 等 ), 无 量 纲 特性 参数 (如 : 浮 升 力 比 N, 
瑞丽 数 Ra, 路 易 斯 数 Le, 格拉 晓 夫 数 Gr 等 ) 对 复合 
方 腔 内 的 自然 对 流 与 热 质 扩散 影响 的 研究 A, 也 
研究 了 复合 方 腔 内 多 孔 介 质 层 与 纯 流体 空间 交界 面 
上 的 应 力 跳 路 系 数 对 方 腔 内 自然 对 流 与 热 质 扩散 的 
影响 站。 对比 研究 了 二 维 与 三 维 复合 方 腔 模型 内 
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自然 对 流 和 热 质 扩散 现象 PON 等 . 以 上 对 自然 对 
流 现 象 的 研究 工作 均 基 于 表征 体 元 (Representative 
Elementary Volume, REV) REE, 采用 有 限 体积 或 者 
有 限 元 等 传统 方法 进行 的 ， 

与 传统 的 模拟 方法 不 同 ， 被 称 为 介 观 方法 的 格 
子 Boltzmann 方法 因 其 计算 效率 较 高 、 编 程 容易 、 边 
界 条 件 处 理 简 单 等 优点 ,被 广泛 应 用 于 微 尺度 流动 
与 换 热 、 多孔 介质 、 晶体 生长 等 一 些 传 统 方 法 难以 胜 
任 的 领域 。 如 张 国 庆 等 [3 采用 格子 Boltzmann 方 
法 在 孔 际 尺度 下 对 放 有 和 多孔 介质 块 的 通道 内 的 流动 
进行 了 模拟 , 研究 了 不 同 多 孔 介 质 模 型 及 雷诺 数 Re 
对 复合 通道 内 流动 的 影响 . Kang 等 03-16] 采用 格子 
Boltzmann 方法 搭建 LB 模型 , 对 孔 际 尺度 下 多 了 筷 介 
质 输 运 过 程 中 的 多 组 分 反应 及 溶质 溶解 过 程 进行 了 
研究 等 . 

由 于 多 孔 介 质 复 杂 物 理 边界 的 无 规则 性 ， 学 者 
的 研究 工作 主要 集中 在 REV 尺度 自然 对 流 的 研究 
上 ,对 孔 际 尺度 下 双 扩 散 混合 对 流 的 研究 较 少 。 因 
此 ， 本 文 利 用 格子 Boltzmann 方法 , 研究 孔 际 尺 度 
PRATT BAL ILD ILE © 分 别 为 0.6, 0.7 和 
0.8 时 , 不 同位 置 多 和 孔 介质 对 顶 盖 驱 动 的 复合 方 腔 内 
的 混合 对 流 与 热 质 扩散 的 影响 ,为 孔 际 尺度 下 多 筷 
介质 内 的 热 质 交 换 的 进一步 研究 提供 一 定 参 考 ， 


图 1 物理 模型 
Fig. 1 Physical Model 


1 物理 与 数学 模型 


本 文 研究 的 物理 模型 如 图 1 Bas. 复合 方 腔 长 
宽 比 H/L=2.0, 方 腔 左 壁 为 高 温 高 浓度 边界 ， 右 壁 
为 低温 低 浓 度 边 界 ， 上 下 壁面 与 多 孔 介 质 固体 颗粒 
边界 均 为 绝热 不 可 渗透 边界 , 顶 盖 以 U=1.0 的 恒定 


陆 RE: 复合 方 腔 顶 盖 驱 动 双 扩散 混合 对 流 格子 Boltzmann 模拟 


速度 水 平 向 右 移动 . 方 腔 内 多 孔 介 质 层 宽 为 D, A 
D/L=1/4, 通过 改变 多 和 孔 介 质 层 与 左 侧 壁 面 间 的 距 
离 Dx, 在 多 和 孔 介质 层 孔 际 率 = 分 别 为 0.6, 0.7 和 0.8 
时 ,探讨 项 盖 驱 动 的 复合 方 腔 内 部 双 扩 散 混 合 对 流 
特性 . 

采用 Wang 等 1718] 提出 的 四 参数 生成 法 (Quar- 
tet Structure Generation Set, QSGS) 构造 孔隙 尺度 
多 和 孔 介质 层 ， 四 参数 分 别 为 生长 核 分 布 率 P., 方向 
生长 率 Py, 相交 互生 长 率 PP 和 固 相 相 含 率 已 , 本 文 
构造 的 多 和 孔 介 质 平均 孔 际 率 e=0.6/0.7/0.8, 平均 固 
体 颗 粒 直 径 Da=18.0, 平均 颗粒 直径 Da = 6/Sv, Sy 
为 多 孔 介 质 固体 颗粒 的 比 表 面积 , Sy = M/A, M 和 
A 分 别 为 二 维 多 孔 介质 固体 颗粒 的 边界 总 周 长 和 面 
HO) ， 所 构造 的 多 孔 介 质 层 如 图 2 所 示 ， 
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图 2 多 和 孔 介 质 层 模型 


Fig. 2 Models of Porous Medium layer 


通过 改变 多 和 孔 介质 层 与 左 侧 壁 面 的 距离 ,研究 
多 了 筷 介 质 层 位 置 对 复合 方 腔 内 热 质 双 扩 散 的 影响 ， 
方 腔 网 格 尺寸 为 800 x 400, 多 和 孔 介质 层 网 格 尺寸 为 
800 x 100， 所 选 5 组 模型 中 多 孔 介 质 层 与 高 温 高 浓 
度 边 界 间 间隔 的 网 格 步 长 如 表 1 所 示 ， 因 本 文 多 筷 
介质 固体 颗粒 边界 为 绝热 不 可 渗透 边界 ， 为 避免 多 
和 孔 介 质 层 对 高 温 高 浓度 与 低温 低 浓 度 边 界 的 有 效 换 
热 面积 的 影响 , 在 模型 D1 与 Ds 中 多 孔 介 质 层 与 左 
右 壁面 间 留 有 5 个 步 长 的 间 际 。 


表 1 多 孔 介质 层 在 方 腔 中 的 位 置 


Table 1 Locations of Porous Medium Layer in 


Cavity 
Dx Dı Do D3 Da Ds 
FBS /lu 5 100 150 200 295 


模型 涉及 的 无 量 纲 特征 参数 : 雷诺 数 Re, 格拉 


nay nj det 
laQA IV Tt 


"cad EE 
f A 
CHN 
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晓 夫 数 Gr, HAR Pr, FEAR Ri, BAM 
数 Le 及 浮 升力 比 NN 定义 如 下 

Re = wt gr Gerth = TLE pp Y 

v’ v2 , a 
Gr Qa Ba(Ch — Ce 

=p ps = TH z 

本 文 模拟 时 普 朗 特 数 Pr 取 定 值 0.7， 格 拉 晓 
夫 数 Gr 取 1.0x104， 路 易 斯 数 Le=2.0， 浮 升力 比 
N=1.0。 复合 方 腔 边界 条 件 设 定 如 下 

EH: U=V=0.0,06=S=1.0 

HE: U =V =0.0,0 = 5S =0.0 

顶板 : U = 1.0, = 0.0,00/0Y = 95/6Y = 0.0 

底板 : U = V =0.0,00/0Y = OS/0Y = 0.0 

多 孔 介 质 男 体 边界 : 00/0n = 0S/0n = 0.0, n 
表示 固体 边界 表面 的 外 法 向 。 


2 热 质 耦 合 双 扩散 Boltzmann 模型 


本 文 在 Guo 等 Bo 提出 的 CLBGK 模型 的 基础 
上 引入 了 深度 分 布 函数 ， 构 造 用 于 模拟 方 腔 内 双 扩 
散 混合 对 流 的 热 质 身 合 的 LBGK 模型 ， 其 密度 分 布 
函数 的 演化 方程 为 
fi(@ + eiði, t + 64) — fi (a, t) = 
— (fi (x,t) — f; (w,t)) /rr + OF; 


AF, e 为 离散 速度 ， 速 度 场 、 温 度 场 及 浓度 场 时 
离散 速度 模型 均 采用 D2Q9 模型 PU, 无 量 纲 松弛 时 
间 re Sie BABE v 满足 关系 uv = h(t — 1/2). F 
衔 态 密度 分 布 函数 fi (wt) 表达 式 为 


Ri 


(2) 


po — 5P/3c? + posi(u) 2 一 0 
Fi" (a,t)=|  P/3c2+posi(u) i=1,2,3,4 (3) 
P/12c? + posi(u) i=5,6,7,8 


silu) = w;[3e;u/c + 4.5(e;u)?/c? — 1.5u?/c?| 
式 中 , MAM wo = 4/9, wi_4 = 1/9, w5_s = 1/36, 
演化 方程 中 外 力 项 F; 表达 式 为 
F; = (1 = 1/27;) wil3 (ce; = u) /© 
+9(e;-u) e;/c?]-F 


假设 方 腔 内 流体 为 不 可 压 牛 顿 流体 ， 忽 略 黏 性 
耗 散 ， 为 便于 处 理由 于 温差 和 浓度 差 而 引起 的 浮 升 
力 项 ， 采 用 Boussinessq 假设 PA, BAD F = 
96r( — To) + g8c(C — Co), 9 为 重力 加 速度 ， 热 
膨胀 系数 Br = —(Op/0T) /oo， 质 膨胀 系数 Bo = 
一 (0p/90) /po. 流 场 宏 观 速度 与 压力 定义 为 


(4) 


8 


8 
u = 5 ceifi P = 区 fit w(t) (5) 
i=1 i=1 


温度 分 布 函数 与 浓度 分 布 函数 的 演化 方程 为 
gi(z 十 eiðs, t + ot) — Gi (x,t) = 


E a 6) 


hi(z 十 eið, t + dt) = hi (x,t) 一 
一 (az 一 (æ, t)) /re 


AF, 无 量 纲 松 儿 时 间 te. To 与 热 扩 散 系数 w、 质 扩 
散 系数 6 满足 关系 


1 1 1 1 
a= J (x 一 5) Cs, B= 3 (x. 一 5) es (8) 


温度 与 浓度 的 平衡 态 分 布 函 数 定义 为 


(7) 


Ser 
g (w, t) =wiT i+ 一 
C 


Ta = ii L 
hod (x,t) = wiC [1 — 3e;- u/c (10) 
宏观 温度 与 宏观 浓度 定义 为 
PS i Bae (11) 
= = 


高 温 高 浓度 壁面 平均 努 塞 尔 数 Nus SPP 
伍德 数 Shay 计算 式 为 


H 
Nisa -5f wi y 而 
HJ, VI lza 
Lf" fae 
ere i pil Y 1 


方 腔 四 周 壁 面 边界 处 理 采 用 Guo 等 P3 提出 的 
具有 二 阶 精度 的 非 平衡 外 推 格式 ， 多 孔 介 质 固体 颗 
粒 边 界 采 用 反弹 格式 2,， 流 场 初始 密度 po = 1.0, 


3 程序 验证 


本 文 多 孔 介 质 固体 颗粒 边界 为 绝热 不 可 渗透 边 
界 ， 与 方 腔 内 纯 流 体 空间 双 扩 散 的 区 别 仅 在 于 边界 
的 不 同 。 为 对 程序 进行 验证 ， 本 文采 用 热 质 耦合 的 
LBGK 模型 模拟 了 二 维 方 腔 内 纯 流体 空间 的 双 扩 
散 混 全 对流， 计算 了 路 易 斯 数 Le=1.0， 浮 升力 比 
N=1.0, 格拉 晓 夫 数 Gr=10“ 时 方 腔 内 浓度 、 温 度 分 
WE, 高 温 高 浓度 壁面 的 平均 努 塞 尔 数 Nuav 及 平均 
会 伍德 数 Shay 的 值 , 并 与 参考 文献 [25] 所 得 结论 进 
行 了 对 比 . 等 温 线 、 流 线 及 等 浓度 线 分 布 对 比如 图 3 
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等 温度 线 医 


等 浓度 线 图 


图 3 本 文 (上 ) 与 文献 [25]( 下 ) 等 温度 线 、 流 线 和 等 浓度 线 分 布 
图 


Fig. 3 Comparison contours of isotherm, streamline and 
isoconcentration between this article (up) and Ref.[25] (down) 


4.8 


4.4 -一 文献 [25] 
= 本 X 
4.0 


Nuw 
Ww 
nN 


图 4 左 壁 壁面 Nuav( 上 ) 和 Shav( 下 ) 随 Ri 的 变化 对 比 图 
(Le=1.0, N=1.0) 


Fig. 4 Comparison diagrams of Nuav (up) and Shav (down) at 
the left wall changes with Ri(Le=1.0, N=1.0) 


所 示 , 高 温 高 浓度 壁面 平均 努 塞 尔 数 Nusv 及 平均 合 
伍德 数 Shaw 的 对 比 结果 如 图 4 所 示 . 从 图 4 中 可 以 
看 出 本 文 结论 与 参考 文献 25] 采用 有 限 体积 法 模拟 
计算 所 得 结论 吻合 较 好 ， 验 证 了 本 文采 用 的 模型 各 
序 的 准确 性 . 


4 计算 结果 与 分 析 


4.1 不 同 结构 参数 下 方 腔 内 热 质 扩散 分 布 特性 

图 5(a) 所 示 为 孔隙 率 e=0.6/0.7/0.8, AIA FL 
介质 层 位 置 下 , 方 腔 内 热 质 扩 散 及 流 线 分 布 图 . 从 图 
5(a) 中 的 流 线 分 布 可 以 看 出 ， 自 由 空间 靠近 顶 盖 位 
置 形成 明显 的 涡流 ， 多 孔 介 质 层 阻力 较 大 ,流体 不 
易 穿 透 。 随 着 多 孔 介 质 层 位 置 右 移 ， 左 侧 自 由 空间 
内 流体 的 自然 对 流 强度 因 浮 升力 作用 逐渐 增强 ,项 
盖 驱 动 引 起 的 局 部 涡流 强度 左 侧 部 分 逐渐 增 大 ,在 
多 孔 介 质 层 右 移 至 位 置 Do 时 ， 涡 流 区 域 相应 演变 
为 上 下 两 个 ;而 方 腔 右 侧 部 分 涡流 逐渐 减 小 ， 在 多 
孔 介 质 层 右 移 至 位 置 Ds 时 中 部 出 现 涡流 ， 由 于 右 
侧 自 由 空间 内 流体 对 流 换 热 强度 随 多 孔 介质 层 的 右 
移 逐 渐 减 弱 , 右 侧 涡流 逐渐 消失 .。 

图 5(b) 所 示 为 孔 际 率 分 别 为 0.7 和 0.8, Ds 位 置 
温度 、 浓度 及 流 线 分 布 图 . 对 比 图 5(a)、(b) 中 的 流 线 
分 布 可 以 发 现 , 孔 际 率 e=0.7 和 e=0.8, 多 孔 介 质 层 
处 于 位 置 Ds 时 , 方 腔 内 流 线 分 布 与 s=0.6( 图 5(a)) 
时 相似 . 且 随 着 孔 际 率 = 的 增 大 , 多 孔 介 质 层 内 渗流 
强度 因 流 动 阻力 的 减 小 而 增强 , FL BRAK ce=0.6 时 位 
置 Ds 中 左 侧 中 部 出 现 的 涡流 与 位 置 Ds 中 右 侧 中 部 
出 现 的 涡流 强度 随 孔 际 率 的 增 大 而 逐渐 增强 ， 且 因 
渗流 强度 的 增强 , 与 另 一 侧 流体 呈 逐 渐 融 合 趋势 ， 

从 图 5(a) 中 的 温度 与 浓度 分 布 去 图 可 以 看 出 ， 
多 孔 介 质 层 对 复合 方 腔 内 的 温度 与 浓度 分 布 影响 显 
著 , 由 于 流动 阻力 的 存在 , 多 孔 介 质 层 内 对 流 强度 比 
自由 空间 弱 很 多 , 很 大 程度 上 限制 了 方 腔 内 的 热 质 
交换 , 使 得 复合 方 腔 内 , 仅 在 左 侧 高 温 高 浓度 与 右 侧 
在 方 腔 中 部 温度 与 浓度 分 布 较为 均匀 。 且 当 多 和 孔 介 
质 层 靠近 高 温 高 浓度 壁面 位 置 (位 置 Di, D2) 时 , 左 
壁 壁面 附近 由 于 对 流 强度 较 弱 ， 传 热 传 质 主要 以 扩 
散 为 主 , 混合 对 流速 度 场 对 温度 与 浓度 分 布 的 影响 
较 小 , 壁面 上 下 温度 与 浓度 梯度 近似 哇 均 匀 分 布 ; 随 
着 多 孔 介 质 层 的 右 移 ， 方 腔 左 侧 混合 对 流 强度 逐渐 
增强 , 对 传 热 传 质 的 影响 逐渐 占据 主导 地 位 , 使 得 左 
壁 壁面 顶板 附近 局 部 温度 与 浓度 分 布 趋 于 均匀 。 对 
比 图 5(a)、(b) 可 以 发 现 , 当 多 孔 介 质 层 位 置 不 变 时 ， 
方 腔 内 混合 对 流 强度 因 和 孔隙 率 的 增 大 而 增强 ， 多 孔 
介质 层 对 方 腔 内 温度 与 浓度 分 布 的 影响 逐渐 减弱 . 
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图 5(a) 多 和 孔 介质 位 置 Dz 不 同时 方 腔 内 温度 (E) 浓度 (下 ) 与 流 线 分 布 图 (e= 0.6) 


Fig. 5(a) Contours of isotherm(up) and isoconcentration(down) and streamline in cavity with different locationsDz (e = 0.6) 
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图 5(b) Ds 位 置 下 孔隙 率 e=0.7 和 e=0.8 时 方 腔 内 温度 (Ze). 浓度 ( 右 ) 与 流 线 分 布 图 


Fig. 5(b) Contours of isotherm(left) and isoconcentration(right) and streamline in cavity with location D3 ate=0.7 and < 一 0.8 
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4.2 不 同 结构 参数 下 高 温 高 浓度 壁面 热 质 交换 特性 

图 6 所 示 为 不 同 孔隙 率 时 ， 顶 盖 驱 动 方 腔 内 多 
孔 介 质 层 不 同位 置 对 左 侧 高 温 高 浓度 壁面 平均 努 塞 
尔 数 Nuay 和 平均 会 伍德 数 Shay 的 影响 曲线 图 .从 
图 中 可 以 看 出 , EME Di~Ds 之 间 ， 当 多 了 筷 介质 层 
和 孔 际 率 = 不 变 时 ， 随 着 多 孔 介 质 层 的 右 移 ， 方 腔 左 
侧 自 由 空间 内 混合 对 流 强 度 逐 渐 增 强 ， 强 化 了 左 壁 
壁面 的 热 质 交换 效率 ， 使 得 方 腔 左 壁 壁面 Nus 和 
Shav 随 多 孔 介质 层 的 右 移 而 增 大 ， 且 和 孔隙 率 © 越 
小 , 增幅 越 明 显 . 在 位 置 Da~Ds 之 间 , 多 孔 介 质 层 
的 存在 对 方 腔 左 侧 混合 对 流 强度 的 影响 减弱 ， 使 得 
多 了 筷 介质 层 的 位 置 Ds KIRE £ 对 Nuav 和 Shav 
的 影响 不 明显 . 仅 当 处 于 位 置 Di t, SANES 
内 的 渗流 强度 因 孔 隙 率 e 的 增 大 而 增强 ， 强 化 了 方 
腔 内 的 热 质 交 换 , 使 得 左 壁 壁面 Nuav 和 Shay 显著 
HK. 
4.3 HERRA Ri 的 影响 

理 查 德 森 数 Ri 是 衡量 混合 对 流 中 强制 对 流 与 
自然 对 流 相对 强 弱 的 关键 因素 ， 为 探究 Ri 对 复 
合 方 腔 壁面 对 流传 热 传 质 特性 的 影响 ， 本 文 对 不 同 


Dı D, D; Dy Ds 


Ri(10-3 <Ri<10*) 下 ， 多 孔 介 质 层 处 于 位 置 Ds 
时 ,复合 方 腔 内 的 对 流传 热 传 质 进行 了 模拟 。 如 图 
7 所 示 为 孔 际 率 ce=0.6 不 同 Ri 下 方 腔 内 的 流 线 分 
ME. 

从 图 中 可 以 看 出 ，Ri 对 复合 方 腔 内 流 场 分 布 的 
影响 显著 ， 当 方 腔 内 为 自然 对 流 (Ri<107t) 时 , 复 
合 方 腔 内 流 线 分 布 较为 复杂 ， 受 顶板 影响 ， 多 孔 介 
质 层 左右 两 侧 顶 板 附近 均 存 在 顺 时 针 方 向 涡流 ， 而 
在 右 侧 中 部 局 部 出 现 道 时 针 方向 涡流 ; 当 流 体 处 于 
混合 对 流 流 态 (107! <Ri<10) 时 , 局 部 涡流 处 于 顶 
板 附 近 ， 逆 时 针 方 向 涡流 消失 ; 当 流 体 处 于 强制 对 
流 流 态 (Ri>10) F, 方 腔 两 侧 流 线 分 布 ， 两 侧 涡 流 
中 心 分 别处 于 方 腔 两 侧 中 部 且 涡 流 占 据 整 个 自由 空 
E. 显然 , 随 着 Ri 的 增 大 , 多 和 孔 介 质 层 内 渗流 强度 
明显 减弱 , 流体 在 多 孔 介 质 层 两 侧 分 别 形成 涡流 , 在 
这 很 大 程度 上 限制 了 方 腔 左右 两 侧 流 体 间 的 热 质 交 
He, 从 而 影响 壁面 对 流传 热 传 质 效率 . 如 图 8 所 示 ， 
不 同和 孔 际 率 下 方 腔 左 壁 壁面 Nus, 和 Shay 随 Ri 变 
化 曲线 图 能 直观 反映 这 一 点 , 方 腔 左 壁 壁面 Nua 和 
Shay BE Ri 的 增 大 逐渐 减 小 . 


图 6 不 同 孔 阶 率 下 Nuav (a) 和 Shay (b) 随 多 孔 介 质 层 位 置 变化 关系 图 


Fig. 6 Nuav (a) and Shay (b) of the left wall changes with porous medium layer locations at different porosity 
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Ri=10-3 i=1. i Ri=10? 
图 7 Ds 位 置 不 同 Ri(10~3 <Ri<103) 下 复合 方 腔 内 流 线 分 布 (==0.6) 
Fig. 7 Streamlines in cavity with different Richardson number (1073 <Ri<10°) at D3 (e=0.6) 
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图 8 位 置 Ds 下 左 壁 壁面 Nuav(a) 和 Shav(b) 随 理 查 德 森 数 
Ri 变化 关系 图 
Fig. 8 Nuav (a) and Shav(b) of the left wall changes with 
Richardson number Ri at different porosity 


5 结 论 


本 文采 用 Boltzmann 方法 在 孔 际 尺 度 下 ， 对 
路 易 斯 数 Ze=2.0， 格 拉 晓 夫 数 Gr=10*, WITH tk 
N=1.0, 普 朗 特 数 Pr=0.7 时 , 多孔 介质 层 与 纯 流 体 
空间 组 成 的 复合 方 腔 内 顶 盖 驱动 的 双 扩 散 混合 对 流 
现象 进行 了 数值 模拟 , 分 析 讨 论 了 多 了 筷 介 质 层 孔 际 
率 、 不 同位 置 及 Ri 对 方 腔 内 热 质 双 扩 散 特 性 的 影 
w, 得 到 结论 如 下 : 

1) 多 和 孔 介质 层 的 存在 对 方 腔 内 项 盖 驱 动 下 的 双 
扩散 特性 影响 显著 ， 且 孔 际 率 越 小 ， 多 和 孔 介 质 层 越 
靠近 高 漫 高 浓度 壁面 , 影响 越 明显 ， 

2) 当 多 和 孔 介质 层 孔 际 率 = 不 变 时 ， 在 位 置 
Di~D3 之 间 ， 方 腔 高 温 高 浓度 壁面 Nua 和 Shay 
随 多 孔 介 质 层 的 右 移 而 增 大 , 在 位 置 D3~Ds 之 间 ， 
多 和 孔 介 质 层 位 置 对 高 温 高 浓度 壁面 Nuav 和 Shay 影 
响 不 显著 .。 


3) 不 同位 置 下 , 孔 际 率 对 高 温 高 浓度 壁面 Nuav 
和 Shay 影响 不 具 统 一 规律 性 . 

4) 当 多 和 孔 介 质 层 处 于 位 置 D3 时 , 高 温 高 浓度 
壁面 Nuav 和 Shay 随 Ri 的 增 大 而 减 小 。 
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